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I.  INTRODUCTION 


An  important  relationship  used  commonly  in  magnetohydrodynamic  (MHD)  descrip¬ 
tions  of  convection  in  the  earth’s  plasma  sheet  is  one  that  links  plasma  pressure,  P,  and 
mass  density,  p.  The  equation  of  state,  P  =  ap^^^,  where  o  is  a  constant  of  motion, 
emerges  from  the  combined  MHD  mass,  momentum,  and  energy  equations  if  a  scalar  pres¬ 
sure  and  an  adiabatic  (constant  entropy)  flow  are  assumed  (see  Siscoe  [1983]).  For  more 
general  conditions,  the  equation  of  state  can  be  used  in  its  polytropic  form,  P  =  op"^, 
where  7  is  the  polytropic  index. 

We  address,  in  this  report,  the  “effective”  polytropic  equation  of  state  in  the  earth’s 
magnetotail  by  developing  further  a  theoretical  model  of  plasma  sheet  convection.  The 
model  extends  the  analysis  of  Kivelson  and  Spence  [1988],  who  showed  that  during  suf¬ 
ficiently  slow  earthward  convection,  a  flux  tube’s  integrity  is  violated  because  cross-tail 
(gradient/curvature)  drift  of  the  energetic  particle  population  is  large,  relative  to  the  bulk 
earthward  (E  xB)  drift.  [This  effect  can  be  incorporated  as  a  cross-tail  heat  flux  in  an 
MHD  description  (Goertz,  personal  communication,  1989).]  Consequently,  pressures  cal¬ 
culated  in  a  magnetotail  of  finite  width  can  be  reduced  significantly  from  those  calculated 
using  two-dimensional  (2D),  lossless  adiabatic  convection  models.  In  this  report,  we  deter¬ 
mine  the  effect  of  the  finite  tail  width  on  number  density  and  find  that  flux-tube  content 
is  not  a  constant  of  motion  along  the  midnight  meridian.  Finally,  we  use  the  modeled 
pressures  and  number  densities  to  define  an  effective  polytropic  equation  of  state. 
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II.  CONVECTION  MODEL 


We  adopt  the  modified  model  of  plasma  sheet  convection  presented  first  by  Erickson 
[1985]  and  discussed  more  recently  by  Kivelson  and  Spence  [1988].  We  assume  for  sim¬ 
plicity  a  zero  dipole  tilt  and  use  geocentric  solar  magnetospheric  (GSM)  coordinates.  At 
a  large  down-tail  distance  Xq  (taken  to  be  -60  Re),  we  assume  that  a  uniform  plasma, 
characterized  by  a  Maxwellian  temperature,  To,  and  a  number  density,  no,  is  the  source  of 
plasma  for  the  tail  plasma  sheet.  The  plasma  is  distributed  uniformly  across  a  tail  of  width 
2R.  In  this  treatment,  we  do  not  include  low  latitude  boundary  layer  (LLBL)  particles 
that  drift  in  from  the  dawn  magnetopause.  The  LLBL  constitutes  a  source  of  additional 
particles  for  the  near-earth  plasma  sheet,  so  by  ignoring  it,  we  underestimate  the  density 
and  pressure  along  the  midnight  meridian.  We  assume  that  earthward  convection  of  the 
plasma  is  driven  by  a  uniform  magnetotail  electric  field,  E,  pointing  from  dawn  to  dusk  and 
that  the  2D  magnetic  field  is  given  by  the  midnight  meridian  cut  through  the  Kp  =  0,  O'*' 
Tsyganenko  [1987]  (T87)  model  (see  his  Figure  2).  Hence  the  model  is  pertinent  to  the 
magnetotail  under  geomagnetically  quiescent  conditions.  Furthermore,  we  assume  that 
pitch  angle  scattering  maintains  isotropic  distributions  upon  convection.  We  describe  the 
motion  of  the  bounce  orbit  at  the  equatorial  (X-Y)  plane  at  Z  =  0  using  bounce-averaged 
drift  theory.  At  a  location  earthward  of  the  source,  (Xe,0,0),  the  spectrum  of  particles 
arriving  from  the  source  cuts  off  above  a  critical  energy  which,  by  conservation  of  energy, 
can  be  related  to  a  critical  energy,  Wo"*,  in  the  source  spectrum.  The  value  of  Wq"'  de¬ 
pends  on  the  parameters  of  the  convection  process  a.3  Wo"*  =  qER/ (^[Vo/V{Xe)]^  —  1^, 

where  V’(Xe)  [Vo]  is  the  specific  volume  of  the  flux  tube  intersecting  Xg  [A^’o]  [Kivelson 
and  Spence,  1988].  Particles  starting  from  (.Yo,V,0)  with  Wo  >  Wq"'  acquire  cross-tail 
velocities  that  sweep  them  duskward  of  midnight  before  they  reach  A’e. 

Macroscopic  plasma  parameters  are  obtained  from  velocity  moments  of  the  distribution 
function,  and  Liouville’s  theorem  specifies  the  mapping  of  phase  sj^ace  density  from  A'o  to 
Xe-  The  upper  bound  to  the  velocity  integrals  that  yield  the  desired  moments  is  defined 
by  Wo"*  [uo"*  =  (2Wo"*/mp)^]  and  is  finite  if  /?  <  oo.  Note  that  we  assume  protons,  with 
mass,  rrip,  dominate  the  plasma  sheet  plasma,  and  we  neglect  electrons  and  heavy  ions. 
Erickson  [1985]  and  Kivelson  and  Spence  [1988]  showed  that  the  plasma  pressure  at  A%  in 
a  finite  width  tail  is  reduced,  relative  to  that  in  an  infinitely  wide  (2D)  tail,  by  the  fraction 

=er/(C)-2C(3  +  2(;^)exp(-C')/3(»)»  (1) 

where  erf(Q  is  the  error  function  and 

C'  =  Wo"*/fcTo  =  r/([Vo/V(Xe)]^  -  l)  (2) 

r  =  qER/kTo  (3) 

Poo  was  obtained  from  the  relation  Poo  oc  V“‘>',  where  the  polytropic  index,  7,  was  taken 
to  be  |.  Kivelson  and  Spence  [1988]  evaluated  Eq.  (1)  and  showed  that  for  low  values  of  r 
(r  <  10),  the  finite  tail  width  effect  can  lead  to  pressure  ratios  significantly  less  than  1  and 
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may  provide  a  rcsohition  to  the  “pressure-balance  inconsistency”  described  by  Erickson 
and  M^o//[1980]. 

We  next  consider  the  degree  to  which  number  density  profiles  arc  modified  by  the  finite 
tail  width  effect.  For  the  case  of  2D,  lossless  convection,  flux-tube  content  is  a  constant  of 
motion,  and  this  requirement  specifies  the  radial  variation  of  n  according  to  nV  =  coyisfant. 
We  next  determine  the  variation  of  n  in  a  finite  width  tail  by  calculating  the  zeroth  velocity 
moment  of  the  distribution  function  whose  second  moment  yielded  Eq.  (1).  Here  we 
assume  that  protons  are  the  dominant  ion  species,  and  we  relate  the  mass  density  to  the 
ion  number  density  by  /?„,  =  nrrip.  We  find  that  the  ratio  of  the  number  density  in  a  tail 
of  finite  width  relative  to  that  in  an  infinitely  wide  tail  (n<3o  =  rioVo/V)  is  given  by 

=  e^/(C)  -  2Ce3-p(-C^)/7r^  (4) 

In  Figure  1,  we  have  plotted  the  number  flen.sity  ratio,  versus  the  flux-tube 

volume  ratio,  Vq/V,  as  a  function  of  r.  As  r  approaches  infinity  (i.e.,  a  2D  magnetotail), 
~  R^oo-  The  departure  of  the  number  density  from  that  predicted  by  constant  flux- 
tube  content  is  less  than  10%  for  r  >  50  (see  tojunost  curve  in  Figure  1)  and  is  only 
weakly  dependent  on  the  compression  ratio,  Vq/V.  On  the  other  hand,  as  r  decreases  to 
small  values,  nn/n^  becomes  significantly  less  than  1  for  sufficiently  large  compression. 
For  r  =  10,  n/j  is  approximately  0.4  (0.6,  0.8)  the  value  j)rcdicted  by  the  requirement  of 
constant  flux-tube  volume  for  a  compression  ratio  of  45  (17,  10). 

To  show  more  clearly  how  flux-tube  contc'iit  change's  with  inward  conve'ction,  w('  have 
plotted  nV  versus  Xe  in  Figure  2a.  The  flux-tid)e  volume  (obtained  from  the  T87  model) 
and  the  number  density  have  been  normalized  by  their  values  at  CO  Ri-j:  Vo  =  7.5  x  lO'® 
m^/Wb  and  Mq  =  0.1  cm~'^  (consistent  with  the  quiet  time  values  observed  near  lunar  orbit 
by  Rich  et  al.  (1973]).  Figure  2a  shows  that  nV  is  nearly  constant  beyond  ~  30  /?£-  for 
all  interesting  values  of  r.  At  these  distances,  the  finite  tail  width  effect  is  unimportant. 
At  radial  distances  less  than  30  /?£,  the  normalized  flux-tube  content  drops  noticeably 
below  1  at  a  location  that  shifts  earthwar<l  with  increasing  r.  For  example,  for  r  <  10, 
nV  =  constant  overestimates  the  density  in  the  near-magnetotail  by  a  factor  >  3. 

In  Figure  2b  we  have  plotted  n/j  versus  distance  along  the  midnight  meridian  as  given 
by  Eq.  (4).  Curves  are  parameterized  by  r  and  arc  normalized  to  0.1  cm“'^  at  GO  Re- 
As  noted  in  Figure  2a,  the  finite  tail  width  has  little  effect  on  the  radial  profile  of  number 
density  beyond  30  Re',  however,  the  dependence  of  n/j  on  r  is  strong  for  A'^  <  20  Re- 
We  know  that  the  quiet  time  central  plasma  sheet  number  density  rarely  exceeds  1  cm“^ 
between  10  and  ~20  Re  [Daximjohann  and  Paschmann,  1989;  Huang  et  ai,  1989].  This 
empirical  constraint,  when  combined  with  the  curves  in  Figure  2b,  argues  for  r  <  7  during 
quiet  times.  We  note  that  r  <  7  is  consistent  with  the  value  determined  by  Kivclson  and 
Spence  [1988],  who  found  that  r  <  5  was  required  for  the  modeled  pressure  profile  in  the 
near-earth  plasma  sheet  not  to  exceed  the  lobe  magnetic  pressiirc.  These  limits  on  r  may 
require  modification  when  the  LLBL  source  of  plasma  is  included  in  the  model. 
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Figure  1.  The  ratio  {njifn^)  of  the  number  density  along  the  midnight  meridian  in  a 
tail  of  width  2R  to  the  number  density  in  a  tail  of  infinite  width  versus  the 
ratio  {Vq/V)  of  the  flux-tube  volumes  at  the  source  location  and  closer  to  the 
earth.  The  curves  are  labeled  t,  the  ratio  of  the  potential  drop  across  half  of 
the  tail  to  the  thermal  energy  per  unit  charge  in  the  source  distribution. 
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l^igure  2.  Normalized  flux-tube  content  (a),  number  d<'nsity  (/<),  iuid  tlu'  “effective”  lo¬ 
cally  determined  polytropic  index  (c)  plotted  isus  diht  iuice  in  the  magnetic 
equatorial  plane  along  the  midnight  meridian  as  a  function  of  r.  In  a,  nV  is 
normalized  by  its  value  at  60  Re-  In  6,  the  em  vi-  is  obtained  by  the  as¬ 
sumption  of  lossless  convection  in  a  tail  of  infinite  width  (i.e.,  7)  V'  =  constant). 
In  c,  for  large  values  of  r  and/or  large  down-tail  distances,  7  is  nearly  for 
r  <  15  and  for  relatively  small  down-tail  distances,  7  may  take  on  values 
much  less  than  |. 
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III.  DETERMINATION  OF  THE  RADIAL  VARIATION  OF  THE  POLYTROPIC  INDEX 


In  this  section,  we  use  the  modeled  pressure  and  number  density  to  define  locally  a 
polytropic  index,  7,  where  P/Pq  =  (n/no)'''.  For  an  infinitely  wide  magnetotail,  r  =  00 
and  7  =  |.  Results  for  other  values  of  r  are  shown  in  Figure  3,  where  we  have  plotted 
versus  Pr/Pq  on  logarithmic  scales.  All  curves  were  normalized  to  a  pressure  and 
density  at  Xq  consistent  with  observations.  For  an  assumed  polytropic  relation  between 
Pr/Pq  and  ur/tio,  the  slopes  of  the  curves  in  Figure  3  determine  local  '^alues  of  7  (dashed 
lines  identify  slopes  corresponding  to  7  =  |  and  1).  The  slope  of  the  topmost  solid  curve 
(r  =  00)  is  |,  corresponding  to  adiabatic  convection  in  the  2D  fluid  limit.  Note  that 
as  T  decreases,  the  slopes  of  the  curves  become  shallower  as  the  density  increases.  This 
indicates  that  flux  tubes  close  to  the  inner  edge  of  the  tail  which  are  greatly  compressed 
have  local  values  of  7  <  |. 

Using  the  slopes  of  the  curves  in  Figure  3,  we  have  obtained  values  of  the  locally  defined 
7  as  a  function  of  distance,  and  we  show  the  results  fin  different  values  of  r  in  Figure  2c. 
Dashed  horizontal  lines  identify  7  =  |,  1,  and  1.  For  A^’e  >  40  Rr,  7  for  all  values 
of  T  that  we  consider.  At  these  large  down-tail  distances,  the  finite  tail  corrections  to  P 
and  n  are  relatively  unimportant,  because  cross-tail  drifts  are  small  compared  to  s\mward 
E  X  B  drifts  and  the  integrity  of  plasma  within  a  flux  tube  is  effectively  maintained  upon 
convection.  However,  for  Xg  <  40  Rr,  7  is  sensitive  to  both  Xe  and  r.  For  small  values 
of  T,  7  departs  from  the  t  —  00  curve.  The  location  where  7  deviates  significantly  from 
I  moves  down  tail  (larger  X^)  with  decreasing  r.  We  note  that  for  nominal  t  (<  7),  our 
model  predicts  7  <  |  at  A'e  ~  10  Re  and  1  <  7  <  |  at  A’e  ~  20  Rr. 
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Figure  3.  Solid  curves  show  the  relationship  between  Pr/Po  and  ur/tiq  on  a  logarithmic 
scale  as  a  function  of  r.  The  curves’  slopes  define  locally  a  polytropic  index,  7; 
the  dashed  lines  indicate  7’s  of  ~  and  1.  As  r  ni)proaches  infinity,  7  approaches 
|;  for  smaller  values  of  r,  slopes  are  shallower  and,  therefore,  7  <  |. 
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IV.  QUALITATIVE  COMPARISON  OF  7  WITH  IN  SITU  DETERMINATIONS 


Recently,  two  indepcrulent  data  surveys  determined  7  empirically  in  the  (piiet-time, 
near-tail  plasma  sheet.  Baumjohann  and  Pa^chmann  [1989]  analyzed  an  extensive  data 
set  compri.sed  of  plasma  pressure  and  density  moments  obtained  at  4.5  s  re.solution  with  the 
ion  release  module  (IRM)  threc'-dimensional  (3D)  plasma  instrument.  Data  Vi/ere  obtained 
when  the  IRM  spacecraft  was  in  the  geomafrnetic  tail  (  —  9  >  XcsM  >  —13,  |lr;5A/|  ^  12, 
2  <  2/vs  ^  —3),  and  individual  measurements  wen*  identified  as  being  either  from  the 
central  plasma  sheet  (CPS)  or  from  the  plasma  sheet  boundary  layer  (PSDL). 

Baumjohann  and  Paschmann  [1983]  plotted  log  pressure  versus  log  density  (both  CPS 
and  PSBL)  in  the  same  format  as  our  Figur<‘  3  and  j)erformed  linear  regression  analysis 
oil  the  data.  They  found  that  the  best-fit  slop*'  implied  7  =  1.C9.  Tlu'y  al.so  invest igat<‘d  7 
in  the  CPS  and  the  PSDL  s('parately  and  found  that,  while  P  and  ri  were  wi-ll  correlated 
and  7  =  1.53  in  the  PSDL,  the  correlation  was  too  low  to  deti'iinine  a  meaningful  index 
in  the  CPS.  Desjute  the  simplicity  of  our  mo<l<'l,  we  believi-  these  results  are  (pialit atively 
consistent  with  its  predictions.  Measurenu'uts  obtained  in  the  PSDL  (i.e.,  at  higher  mag¬ 
netic  latitudes,  where  the  vertical  distance  from  the  neutral  sheet  is  relatively  large)  will 
normally  map  to  ecpiatorial  locations  well  down  tail  from  IRM’s  XasM  location.  Finite 
tail  width  effects  are  relatively  small  at  the.se  distances  and,  th('r<'for(',  we  predict  that 
7  will  differ  little  from  On  the  other  han<l,  nu'asurenu'iits  obtaiiu'd  in  the  CPS  are 
presumably  taken  much  closer  to  the  magnetic  e([uator  and  theri'fore  ma])  to  eiiuatorial 
locations  in  the  near-vicinity  of  the  IRM  orbit.  At  these  distances,  our  model  predicts  7’s 
ranging  from  1  to  |,  with  values  highly  d<'peudent  on  both  location  and  t.  Consequently, 
a  poor  correlation  between  CPS  P  and  n  is  expected  in  a  statistical  sample  that  has  not 
been  binned  according  to  location  and  geomagnetic  activity.  Finally,  Baumjohann  and 
Paschmann  [1383]  binned  the  combined  CPS  and  PSBL  data  into  two  auroral  electrojet 
(AE)  index  bins  (  ^  100  nT)  and  found  7  ~  f  for  the  active  subset  and  ~  |  for  the 
geomagneticcilly  quieter  subset.  If  larger  E  and  r  correspond  to  higher  levels  of  AE,  these 
results  are  qualitatively  consistent  with  the  model  predictions  and  are  also  quantitatively 
plausible.  Binning  of  the  data  by  location  and  activity  would  jirovide  additional  tests  of 
this  model  and  forthcoming  modifications  of  it  that  will  include  LLBL  plasma  sources. 

Huang  e.t  al.  [1389]  also  determined  7  empirically  in  tlii'  magnetotail  plasma  sheet.  Their 
survey  is  comprised  of  plasma  density  and  temperature  {kT  =  P/it)  estimates  obtained  by 
the  low  energy  proton  and  eh-ctron  differential  energy  analyzers  (LEPEDEA)  instrument 
on  the  international  sun-earth  explorer- 1  (ISEE-1)  siiacecraft.  The  Huang  et  al.  study 
concentrated  on  measurements  from  geomagnetically  (juiet  intervals  [AE  <100  iiT)  when 
ISEE-1  was  in  the  CPS  (between  —10  <  XasM  ^  —23  /?/•;).  Log  temperature  was  ])lotted 
versus  log  number  density,  and  the  slope  of  th<*  best-fit  linear  trend  was  used  to  determine 
7.  Even  though  the  correlation  coefficient  between  temperature  and  density  is  n'latively 
low,  Huang  et  al.  concluded  that  the  two  quantities  are  anticorrelated  and  thendore  that 
7  <  1.  In  our  model,  we  find  that  7  can  be  signifit  antly  <  |  in  the  slowly  convecting  near¬ 
tail  CPS  but  that  it  is  everywhere  >  1.  Elsewhere  [Spr.nr.r.  and  Kivr.lson,  in  preparation, 
1990),  we  will  show  that  even  for  a  locally  defined  7  >  1,  the  temperature  can  increase 
with  distance  in  regions  where  the  density  is  decreasing  with  distance. 
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V.  SUMMARY 


We  have  developed  a  theoretical  model  of  plasma  sheet  convection  for  a  plasma  source 
in  the  distant  tail  to  determine  how  number  density  varies  with  down-tail  distance  near 
the  midnight  meridian.  We  find  that  fiux-tube  content  is  not  conserved  during  periods  of 
slow  convection.  This  means  that  an  MHD  formulation  of  the  convection  problem  may 
have  to  include  not  only  the  previously  noted  cross-tail  heat  flux  but  also  a  cross-tail 
mass  flux.  Modeled  number  densities  were  combined  with  previously  modeled  pressures 
to  construct  locally  a  polytropic  index.  We  find  that  7  ranges  between  1  and  |.  Low 
values  of  7  occur  during  periods  of  small  r  (intervals  of  slow  convection)  in  the  near-tail 
plasma  sheet  (<  20  Re)i  while  values  nearer  |  occur  during  periods  of  larger  r  (intervals 
of  significantly  enhanced  convection)  and/or  at  larger  down-tail  distances  (>  30  Re)-  We 
feel  that  this  model  may  in  part  explain  some  of  the  apparent  discrepancy  in  the  two 
empirical  determinations  [Baumjohann  and  Paschmann,  1989;  Huang  et  a/.,  1989]  of  the 
plasma  sheet  7.  Further  meaningful  comparison  with  data  may  require  refinement  of  the 
model. 

Possibly  the  most  serious  flaw  in  the  present  work  is  the  neglect  of  an  LLBL  source. 
The  contribution  of  such  a  source  will  be  described  elsewhere.  Preliminary  calculations 
show  that  the  LLBL  source  modifies  the  solution  along  the  midnight  meridian  by  increas¬ 
ing  n  and  P  in  regions  where  the  pres.sure  from  the  distant  tail  source  drops  below  the 
lobe  magnetic  pressure.  Including  the  LLBL  source  should  reduce  the  pressure  imbalance 
apparent  in  Figure  2  of  Kivthon  and  Spence  [1988]  for  small  r  inside  of  20  Re-  The  LLBL 
source  is  also  capable  of  partially  poptilating  the  high-energy  part  of  the  distribution  in 
the  near-earth  tail  that  is  absent  (becau.se  of  the  energy  cutoff  previously  described)  if  only 
a  tail  source  is  considered. 

Features  of  the  model  that  require  Birther  examination  include  the  variation  of  the 
moments  with  cross-tail  distance  and  the  effects  of  skewed  or  nonviniform  electric  fields 
lAtkinson,  1984]  and  of  3D  magnetic  structiire.  The  consideration  of  y-dependence  raises 
the  question  of  electron  reponse.  The  model  implies  a  y-dependence  of  Tf/Ti  in  the  near- 
earth  pleisma  sheet,  but  no  striking  variation  has  been  reported.  Possibly  the  consideration 
of  the  field-aligned  currents  consistent  with  this  model  will  eliminate  this  apparent  discrep¬ 
ancy  between  model  and  data. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  fiinctions  as  an  “architect-engineer”  for  national  security 
projects,  specializing  in  advanced  military  space  systems.  Providing  research  support,  the 
corporation’s  Laboratory  Operations  conducts  experimental  and  theoretical  investigations  that 
focus  on  the  ai^lication  of  scientific  and  technical  advances  to  such  systems.  Vital  to  the  success 
of  these  investigations  is  the  technical  staffs  wide-ranging  expertise  and  its  ability  to  stay  current 
with  new  developments.  This  expertise  is  enhanced  by  a  research  program  aimed  at  dealing  with 
the  many  problems  associated  with  rapidly  evolving  space  systems.  Contributing  their  capabilities 
to  the  research  effort  are  these  individual  laboratories: 

Acrophysics  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat  transfer 
and  flight  dynamics;  chemical  and  electric  propulsion,  propellant  chemistry,  chemical 
dynamics,  environmental  chemistry,  trace  detection;  spacecraft  structural  mechanics, 
contamination,  thermal  and  structural  control;  high  temperature  thermomechanks,  gas 
kinetics  and  radiation;  cw  and  pulsed  chemical  and  excipier  laser  development, 
including  chemical  kinetics,  spectroscopy,  optical  resonators,  beam  control,  atmos¬ 
pheric  propagation,  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions,  atmospheric 
optics,  light  scattering,  state-spedflc  chemical  reactions  and  radiative  signatures  of 
missile  plumes,  sensor  out-of-feM-of-view  rejection,  applied  laser  spectroscopy,  laser 
chemistry,  laser  optoelectronics,  solar  cell  physics,  battery  electrochemistry,  space 
vacuum  and  radiation  effects  on  materials,  lubrication  and  surface  phenomena, 
thermionic  emission,  photosensitive  materials  and  detectors,  atomic  frequency  stand¬ 
ards,  and  environmental  chemistry. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device  physics, 
compound  semiconductors,  radiation  hardening;  electro-optics,  quantum  electronics, 
solid-state  lasers,  optical  propagation  and  communfeations;  microwave  semiconductor 
devices,  microwave/millimeter  wave  measurements,  diagnostics  and  radiometry,  micro¬ 
wave/millimeter  wave  thermionic  devices;  atomic  time  and  frequency  standards; 
antennas,  rf  systems,  electromagnetic  propagation  phenomena,  space  communication 
systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals,  alloys, 
ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  nondestructive 
evaluation,  component  failure  analysis  and  reliability;  fracture  mechanics  and  stress 
corrosion;  analysis  and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures 
as  well  as  in  space  and  enemy-induced  environments. 

Space  Sciences  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray  physics, 
wave-partkle  interactions,  magnetospheric  plasma  waves;  atmospheric  and  ionospheric 
physics,  density  and  composition  of  the  upper  atmosphere,  remote  sensing  using 
atmospheric  radiation;  solar  physics,  infrared  astronomy,  infrared  signature  analysis; 
effects  of  solar  activity,  magnetic  storms  and  nuclear  explosions  on  the  earth’s 
atmosphere,  ionosphere  and  magnetosphere;  effects  of  electromagnetic  and  particulate 
radiations  on  space  systems;  space  instrumentation. 


